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Abstract 

We review analytic methods to perform the post-Newtonian expansion of 
gravitational waves induced by a particle orbiting a massive compact body, 
based on the black hole perturbation theory. There exist two different meth- 
ods of the post-Newtonian expansion. Both are based the Teukolsky equation. 
In one method, the Teukolsky equation is transformed into a Regge- Wheeler 
type equation that reduces to the standard Klein-Gordon equation in the flat 
space limit, while in the other method, which were introduced by Mano, Suzuki 
and Takasugi relatively recently, the Teukolsky equation is directly used in its 
original form. The former has an advantage that it is intuitively easy to un- 
derstand how various curved space effects come into play. However, it becomes 
increasingly complicated when one goes on to higher and higher post-Newtonian 
orders. In contrast, the latter has an advantage that a systematic calculation 
to higher post-Newtonian orders is relatively easily implementable, but other- 
wise so mathematical that it is hard to understand the interplay of higher order 
terms. In this paper, we review both methods so that their pros and cons may 
be clearly seen. We also review some results of calculations of gravitational 
radiation emitted by a particle orbiting a black hole. 



1 Introduction 
1.1 General 

In the past several years, there has been substantial progress in the projects of 
ground-based laser interferometric gravitational wave detectors which include 
LIGO [1], VIRGO [2], GEO600 [3], and TAMA300 [75, 9, 4]. TAMA300 has 
been in operation since 1999. LIGO and GEO600 started to operate in 2002. 
VIRGO is expected to start operation in 2003. As another type of detectors, 
resonant bar detectors [8] have been operating continuously for the past several 
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years. There arc several future projects as well. Most importantly, the Laser 
Interferometer Space Antenna (LISA) project is in progress [5, 6]. For a review 
of ground and space laser interferometers, see, e.g., [62]. 

The detection of gravitational waves will be done by extracting gravitational 
wave signals from noisy data stream. In developing the data analysis strategy, 
the detailed knowledge of the gravitational wave forms will help us greatly to 
detect a signal, and to extract the physical information of its source. Thus, it 
has become a very important problem for theorists to predict with sufficiently 
good accuracy the wave forms from possible gravitational wave sources. 

Gravitational waves arc generated by dynamical astrophysical events and 
they are expected to be strong enough to be detected when compact stars such 
as neutron stars or black holes are involved in such events. In particular, co- 
alescing compact binaries arc considered to be the most promising sources of 
gravitational radiation that can be detected by the ground-based laser interfer- 
ometers. The last inspiral phase of a coalescing compact binary, in which the 
binary stars orbit each other for ^ 10^ cycles, will be in the bandwidth of the 
interferometers, and it may be not only detected but also provides us important 
astrophysical information of the system by comparing the inspiral waveform 
with theoretical templates if they have sufficient accuracy. Such an event may 
be detected by the initial phase interferometers if we are lucky, and it will be 
detected by the various advanced detectors [49, 55]. 

To predict the wave forms, a conventional approach is to formulate the Ein- 
stein equations with respect to the flat Minkowski background and apply the 
post-Newtonian expansion to the resulting equations (see the subsection below). 

In this paper, however, we review a different approach, namely the black hole 
perturbation approach. In this approach, binaries are assumed to consist of a 
massive black hole and a small compact star which is taken to be a point particle. 
Hence, its applicability is constrained to the case of binaries with large mass 
ratio. Nevertheless, there are several advantages that cannot be overlooked. 

Most importantly, the black hole perturbation equations take full account of 
general relativistic effects of the background spacetime and they are applicable 
to arbitrary orbits of a small mass star. In particular, if a numerical approach 
is taken, gravitational waves from highly relativistic orbits can be calculated. 
Then, if we can develop a method to calculate gravitational waves to a siiffi- 
ciently high PN order analytically, it can not only give insight into how and 
when general relativistic effects become important, by comparing with numer- 
ical results, but also give us a knowledge complementary to the conventional 
post-Newtonian approach, about such as yet-unknown higher PN order terms 
or general relativistic spin effects. 

Moreover, one of the main targets of LISA is to observe phenomena asso- 
ciated with formation and evolution of super-massive black holes in galactic 
centers. In particular, a gravitational wave event of a compact star spiraling 
into such a super-massive black hole is indeed the case for the black hole per- 
turbation theory. 
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1.2 Post-Newtonian expansion of gravitational waves 

The post-Newtonian expansion of general relativity assumes that the internal 
gravity of a source is small so that the deviation from the Minkowski metric 
is small, and that velocities associated with the source are small compared 
to the speed of light c. When we consider the orbital motion of a compact 
binary system, these two conditions become essentially equivalent to each other. 
Although, both conditions may be violated inside each of the compact objects, 
this is not regarded as a serious problem of the post-Newtonian expansion as 
long as we are concerned with gravitational waves generated from the orbital 
motion, and the two bodies are usually assumed to be point-like objects in the 
calculation. In fact, recently Itoh, Futamasc, and Asada [37, 38] developed a 
new post-Newtonian method which can deal with a binary system in which 
the constituent bodies may have strong internal gravity, based on earlier work 
by Futamasc and Schutz [31, 32]. They derived the equations of motion to 
2.5PN order and obtained a complete agreement with the Damour-Deruelle 
equations of motion [24, 25], which assumes the validity of the point-particle 
approximation. 

There are two existing approaches of the post- Newtonian expansion to calcu- 
late gravitational waves: one developed by Blanchet, Damour, and Iyer [15, 12] 
and another by Will and Wiseman [84] based on previous work by Epstein, 
Wagoner, and Will [27, 81]. In both approaches, the gravitational waveforms 
and luminosity are expanded in time derivatives of radiative multipoles, which 
arc then related to some source multipoles (the relation between them contains 
the "tails") which are expressed as integrals over the matter source and the 
gravitational field. The source multipoles are combined with the equations of 
motion to obtain explicit expressions in terms of the source masses, positions, 
and velocities. 

One issue of the post-Newtonian calculation arises from the fact that the 
post-Newtonian expansion can be applied only to the near zone field of the 
source. In the conventional post-Newtonian formalism, the harmonic coordi- 
nates are used to write down the Einstein equations. If we define the deviation 
from the Minkowski metric as 

/^''^ = ,/^5''^-^?''^ (1) 

the Einstein equations are schematically written in the form 

nh^"" = 16n\g\T^"' + A'*''(/i), (2) 

together with the harmonic-gauge condition, d^h^^" = 0, where □ = rjf^^d^d^, 
is the D'Alambertian operator in flat space time, rj^'^ = diag(— 1, 1, 1, 1), and 
K^'^(h) represents the non- linear terms in the Einstein equations. The Einstein 
equations (2) are integrated using the flat space retarded integrals. In order 
to perform the post-Newtonian expansion, if we naively expand the retarded 
integrals in power of 1/c, there appears divergent integrals. This is a technical 
problem which arises due to the near zone nature of the post- Newtonian approx- 
imation. In the BDI approach, in order to integrate the retarded integrals, and 
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to evaluate the radiative multipole moments at infinity, two kind of approxima- 
tion method are introduced. One is the multipolar post-Minkowski expansion, 
which can be applied to a region outside the source including infinity, and the 
other is the near zone, post-Newtonian expansion. These two expansions are 
matched in the intermediate region where both expansions are valid, and the 
radiative multipole moments are evaluated at infinity. In the WW approach, 
the retarded integrals arc evaluated directly, without expanding in terms of 1/c, 
in the region outside the source in a novel way. 

The lowest order of the gravitational waves are given by the Newtonian 
quadrupolc formula. It is standard to call the post-Newtonian formulas for the 
wave forms and luminosity which contain terms up to 0{{v/c)") beyond the 
Newtonian quadrupole formula as the n/2PN formulas. Evaluation of gravita- 
tional waves emitted to infinity from a compact binary system has been suc- 
cessfully carried out to the 3.5 post-Newtonian (PN) order beyond the lowest 
Newtonian quadrupole formula (apart from an undetermined coefEcient that 
appears at 3PN order) in the BDI approach [15, 12]. The computation of the 
3.5PN flux requires the 3.5PN equations of motion. See a review by Blanchet [13] 
for details on post-Newtonian approaches. Up to now, both approaches give the 
same answer for the gravitational wave forms and luminosity to 2PN order. 

1.3 Linear perturbation theory of black hole 

In the black hole perturbation approach, we deal with gravitational waves from 
a particle of mass /x orbiting a black hole of mass M, assuming /i <SC M. The 
perturbation of a black hole spacetime is evaluated to linear order in /x/M. The 
equations are essentially in the form of Eq. (2) with 77^^ replaced by the back- 
ground black hole metric g^^ and the higher order terms h-{h)^,y neglected. 
Thus, apart from the assumption ji ^ M, the black hole perturbation ap- 
proach is not restricted to slow-motion sources, nor to small derivations from 
the Minkowski spacetime, and the Green function used to integrate the Ein- 
stein equations contains the whole curved spacetime effect of the background 
geometry. 

The black hole perturbation theory was originally developed as a metric per- 
turbation theory. For non-rotating (Schwarzschild) black holes, a single master 
equation for the metric perturbation was derived by Regge- Wheeler [61] for the 
so-called odd parity part, and later by Zerilli [86] for the even parity part. These 
equations have a nice property that they reduce to the standard Klein- Gordon 
wave equation in the flat space limit. However, no such equation has been found 
in the case of a Kerr black hole so far. 

Then, based on the Newmann- Penrose null-tetrad formalism, in which the 
tetrad components of the curvature tensor are the fundamental variables, a 
master equation for the curvature perturbation was first developed by Bardeen 
and Press [11] for a Schwarzschild black hole without source {T^" = 0), and by 
Teukolsky [78] for a Kerr black hole with source {T^^^ ^ 0). The master equation 
is called the Teukolsky equation, and it is a wave equation for a null-tetrad 
component of the Weyl tensor V'o or xjji. In the source- free case, Chrzanowski [21] 
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and Wald [82] developed a method to construct the metric perturbation from 
the curvature perturbation. 

The Teukolsky equation has, however, a rather complicated structure as a 
wave equation. Even in the flat space limit, it docs not reduce to the standard 
Klein-Gordon form. Then, Chandrasekhar showed that the Teukolsky equation 
can be transformed to the form of the Regge- Wheeler or Zerilli equation for 
the source- free Schwarzschild case [19]. A generalization of this to the Kerr 
case with source was done by Sasaki and Nakamura [66, 67]. They gave a 
transformation that brings the Teukolsky equation to a Regge- Wheeler type 
equation that reduces to the Regge- Wheeler equation in the Schwarzschild limit. 
It may be noted that the Sasaki-Nakamura equation contains an imaginary part, 
suggesting that either it is unrelated to a (yet-to-be-found) master equation for 
the metric perturbation for the Kerr geometry or it implies the non-existence 
of such a master equation. 

As mentioned above, an important difference between the black hole pertur- 
bation approach and the conventional post-Newtonian approach appears in the 
structure of the Green function used to integrate the wave equations. In the 
black hole perturbation, the Green function takes account of the curved space- 
time effect on the wave propagation, which implies complexity of its structure in 
contrast to the flat space Green function. Thus, since the system is linear in the 
black hole perturbation approach, the most non-trivial task is the construction 
of the Green function. 

There are many papers that deal with a numerical evaluation of the Green 
function and calculations of gravitational waves induced by a particle. See 
Breuer [18], Chandrasekhar [20], and Nakamura, Oohara and Kojima [48] for 
reviews and for references on earlier papers. 

Here, we are interested in an analytical evaluation of the Green function. 
One way is to adopt the post- Minkowski expansion assuming GM/c? ^ r. 
Note that, for bound orbits, the condition GMjf? <C r is equivalent to the the 
condition for the post-Newtonian expansion jf? <^ \. If we can calculate 
the Green function to a sufficiently high order in this expansion, we may be 
able to obtain a rather accurate approximation of it that can be applicable to 
a relativistic orbit fairly close to the horizon, possibly to the radius as small as 
the inner-most stable circular orbit (ISCO), which is given by rjsco = QGM/(? 
in the case of a Schwarzschild black hole. 

It turns out that this is indeed possible. Though there arise some compli- 
cations as one goes to higher PN orders, they are relatively easy to handle as 
compared to situations one encounters in the conventional post-Newtonian ap- 
proaches. Thus, very interesting relativistic effects such as tails of gravitational 
waves can be investigated easily. Further, we can also easily investigate conver- 
gence properties of the post-Newtonian expansion by comparing a numerically 
calculated exact result with the corresponding analytic but approximate result. 
In this sense, the analytic black hole perturbation approach can provide an 
important test of the post-Newtonian expansion. 
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1.4 Brief historical notes 



Let us briefly review some of past works on post-Newtonian calculations in 
the black hole perturbation theory. Although the literature on numerical cal- 
culations of gravitational waves emitted by a particle orbiting a black hole is 
abundant, there are not so many papers that deal with the post- Newtonian 
expansion of gravitational waves, mainly because it had not been necessary 
until recently when the construction of accurate theoretical templates for the 
interfcrometric gravitational wave detectors had become an urgent issue. 

In the case of orbits in the Schwarzschild background, one of the earliest 
works was done by Gal'tsov, Matiukhin and Petukhov [34], who considered a 
case when a particle is in a slightly eccentric orbit around a Schwarzschild black 
hole, and calculated the gravitational waves up to IPN order. Poisson [56] 
considered a circular orbit around a Schwarzschild black hole and calculated 
the wave forms and luminosity to 1.5PN order at which the tail effect appears. 
Cutler, Finn, Poisson and Sussman [23] also worked on the same problem nu- 
merically by using the least square fitting technique to the numerically evaluated 
data for the luminosity, and obtained a post-Newtonian formula for the lumi- 
nosity to 2.5PN order. Subsequently, a highly accurate numerical calculation 
was carried out by Tagoshi and Nakamura [72]. They obtained the formulas 
for the luminosity to 4PN order numerically by using the least square fitting 
method. They found the \ogv terms in the luminosity formula at 3PN and 4PN 
orders. They concluded that, although the convergence of the post-Newtonian 
expansion is slow, the luminosity formula accurate to 3.5PN order will be good 
enough as theoretical templates for the ground-based interferometers, to rep- 
resent the orbital phase evolution of coalescing compact binaries. After that, 
Sasaki [65] found an analytic method and obtained formiilas which arc needed to 
calculate the gravitational waves to 4PN order. Then, Tagoshi and Sasaki [73] 
obtained the gravitational wave forms and luminosity to 4PN order analytically, 
and confirmed the results of Tagoshi and Nakamura. These calculations were 
extended to 5.5PN order by Tanaka, Tagoshi, and Sasaki [77]. 

In the case of orbits around a Kerr black hole, Poisson calculated the 1.5PN 
order corrections to the wave forms and luminosity due to the rotation of the 
black hole and showed that the result agrees with the standard post-Newtonian 
effect due to spin-orbit coupling [57]. Then, Shibata, Sasaki, Tagoshi and 
Tanaka [69] calculated the luminosity to 2.5PN order. They calculated the lumi- 
nosity from a particle in circular orbit with small inclination from the equatorial 
plane. They used the Sasaki-Nakamura equation as well as the Teukolsky equa- 
tion. This analysis was extended to 4PN order by Tagoshi, Shibata, Tanaka 
and Sasaki [74] in which the orbit of the test particle was restricted to circular 
ones on the equatorial plane. The analysis in the case of slightly eccentric orbit 
on the equatorial plane was also done by Tagoshi [70] to 2.5PN order. 

Tanaka, Mino, Sasaki and Shibata [76] considered the case when a spinning 
particle is in a circular orbit near the equatorial plane of a Kerr black hole, 
based on the Papaetrou equations of motion for a spinning particle [52] and the 
energy momentum tensor of a spinning particle by Dixon [26]. They derived 
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the luminosity formula to 2.5PN order which includes the linear order effect of 
the particle's spin. 

The absorption of gravitational waves into the black hole horizon, appearing 

at 4PN order in the Schwarzschild case, was calculated by Poisson and Sasaki 
for a particle in a circular orbit [59]. The black hole absorption in the case of 
rotating black hole appears at 2.5PN order [33]. Using a new analytic method 
to solve the homogeneous Teukolsky cqiiation found by Mano, Suzuki and Taka- 
sugi [43] , the black hole absorption in the Kerr case was calculated by Tagoshi, 
Mano, and Takasugi [71] to 6.5PN order beyond the quadrupole formula. 

If gravity is not described by the Einstein theory but by the Brans-Dicke 
theory, there will appear scalar type gravitational waves as well as transverse- 
traceless gravitational waves. Such scalar type gravitational waves were calcu- 
lated to 2.5PN order by Ohashi, Tagoshi and Sasaki [51] in the case when a 
compact star is in a circular orbit on the equatorial plane around a Kerr black 
hole. 

In the rest of the paper, we use the units c = G = 1. 



2 Basic formulas for the black hole perturbation 
2.1 Teukolsky formalism 

In terms of the conventional Boyer-Lindquist coordinates, the metric of a Kerr 
black hole is expressed as 



S 
A 



j:de^ 



(3) 



where E = r^+a^ cos^ 6 and A = r^ — 2Mr+o^. In the Teukolsky formalism [78], 
the gravitational perturbations of a Kerr black hole are described by a Newman- 
Penrose quantity ip4 = —Ca0-yS'n"fh'^n'^fh^, where Ca^-ys is the Weyl tensor and 

n" = {{r^ + a%-A,0,a)/i2E), 

m" = (iasin^,0, l,Vsin6l)/(\/2(r + iacos6')). 

The perturbation equation for (j) = p~'^tp4, p = {r — ia cos 0)~^, is given by 

sO(j) = 4n'£f. (4) 

Here, the operator gO is given by 



r(r2 + a2) 



sO = - 



2^2 



a 
A 



+A-''a^(A"+ia^) + -—dg{smede) + 2s 



-2s 



sin 6 



— r — la cos ( 



a(r — M) i cos 



sin 



dt-s{scot^0-l), 



(5) 
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with s = —2. The source term f is given by 
f = 2{B'^ + B*'), 

B'2 = -lp'pL_,[p-'Lo{p-'p-'T^n)] 

^ p'pA^L_,[p-Yj+ip-^p-'A-'T^ 



2V2 



^ p«pAV+[p-VA-iL_i(p-V'7kn)], (6) 



2V2 



where 



TTl 

Ls = do + — — r — awsinO + s cot 6, (7) 
sin^ 

J+ = dr + iK/A; K = {r"^ + 0^)0; - ma , 

and Tnn, and are the tetrad components of the energy momentum 

tensor (T„„ = T^yn^n"^ etc.). The bar denotes the complex conjugation. 

If we set s = 2 in Eq. (4), with appropriate change of the source term, it 
becomes the perturbation equation for tjjo. Moreover, it describes the perturba- 
tion for scalar field (s = 0), neutrino field (|s| = 1/2), and electromagnetic field 
(|s| = 1) as well. 

We decompose V4 into the Fourier-harmonic components according to 

p-V4 = E / du;e-''^'+''^'' -2Sem{0)RimUr), (8) 
em 

The radial function Remu and the angular function sStm{&) satisfy the Teukolsky 
equations with s = — 2 as 

A4fi%^Vnr-)i?w=T,_, (9) 



dr \A dr 
1 d ( ■ n'^ \ 22-2/1 (m- 2cos6') 



sm dO I dti ) sm 



-|- 4aa; cos 9 — 2-\- 2mauj + A 
The potential V{r) is given by 



.^Sem = 0. (10) 



K'^ + 4i(r-M)K , , 

V{r)= ^ + 8ia;r + A, (11) 

where A is the eigenvalue of _2<S'^^. The angular function sSem{9) is called the 
spin-weighted spheroidal harmonic, which is usually normalized as 



p7\ 

Jo 



.2Sem\'^ sm0d0 = 1. (12) 
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In the Schwarzschild limit, it reduces to the spin- weighted spherical harmonic 
with A £{£+ 1). In the Kerr case, however, no analytic formula for A is known. 
The source term T^mw is given by 

/ qau> 
dndtp-'p-\B'^ + s-)e-»-^+-*^H^, (13) 



We mention that for orbits of our interest, which are bounded, T^mw has support 
only in a compact range of r. 

We solve the radial Teukolsky equation by using the Green function method. 
For this purpose, we define two kinds of homogeneous solutions of the radial 
Teukolsky equation: 

r ntrans A 2 -ifer* for r ^ r 

^""^ I C'lmS'^^e''^'' forr^+oo ^ ' 

where k = u — ma/2Mr+ and r* is the tortoise coordinate defined by 



fdr* 



dr 



2Mr+ , r-r+ 2Mr_ r - r_ , , 



where r± = M zL \J M"^ — , and wc have fixed the integration constant. 

Combining with the Fourier mode e^"^*, we see that R^ll^^ has no outcoming 
wave from past horizon, while i?"P has no incoming wave at past infinity. Since 
these are the property of waves causally generated by a source, a solution of 
the Teukolsky equation which has purely outgoing property at infinity and has 
purely ingoing property at the horizon is given by 

(17) 

where the Wronskian Wimu is given by 

We^^ = 2iu:C\ZZ'BT-^. (18) 
Then, the asymptotic behavior at the horizon is 

Dtrans A 2 ^ — iA;r* poo 

(19) 

while the asymptotic behavior at infinity is 



RimUr^<^)^T;-^ dr' """\,;jr^ ^ ^ZZye^-^- . (20) 
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Wc note that the homogeneous Teukolsky equation is invariant under the 
complex conjugation followed by the transformation m —m and w —^ —oj. 
Thus, we can set Ft^^^ = where the bar denote the complex conjuga- 

tion. 

We consider T^^ of a monopole particle of mass /x. The energy momentum 
tensor takes the form, 



S{r-r{t))S{0-e{t))S{^-^{t)). (21) 



T, sin 6dt/dT dr dr 

where = (t, r{t),9{t),ip{t)) is a geodesic trajectory and r = T{t) is the proper 
time along the geodesic. The geodesic equations in the Kerr geometry are given 
by 



dr 
^d^ 



dt 

dr 
^d^ 



C - cos^ ela'^il - E^) + = 0(9), 

^ sin 6* J J 

„2 I „2 



(aE - ^)asm'' e + '^^^(E{r^ + a^) - ah) =T, 
V sin 0/ A \ / 



(22) 



where E, 1^ and C are the energy, the z-component of the angular momentum 
and the Carter constant of a test particle, respectively,^ and 



R = [^(r^ + a^) - ahf - l^[{Ea - hf + + C]. 



(23) 



Using Eq. (22), the tetrad components of the energy momentum tensor are 
expressed as 



Tun = t^^6{r-r{t))S{9-9{t))6{^-^{t)), 
sm9 

Tmn = ti^S{r-r{t))S{9-9{t))6{<fi-ip{t)), 



C'rr 



sm( 



'-6{r-r{t))d{e-9{t))6{^-v{t)), 



(24) 



where 



n 



E{r^+a^)-ak + i:^ 
dr 



^ These constants of motion arc those measured in units of fi. That is, if expressed in 
the standard units, E, Iz and C in Eq. (22) are to be replaced with E/n, Izjp. and C/fj?, 
respectively. 
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Cm n 
Crnrn 



2^/2I]2^ 



CIT 



i sin 6* I aE - 



sm 



d0' 



P 

2Ei 



ismeiaE ^4rA +S 

V sin^ 6// 



d£ 

sin" QJ ' dr 



(25) 



and i = dt/dr. Substituting Eq. (6) into Eq. (13) and performing integration 
by part, we obtain 



Ttr„^ = 4^ / dt [ rf6ie''^*-'™¥'(*) 
\/27r J-oo J 



--L\{p-*Ll{p^S)}Cnnp-^p-'6{r - r{t))5{9 - 0{t)) 



+ ^^{LlS + ia{p - p) sin 6S) J+ \ -^S{r - r{t))S{9 - 9{t)) 



V2p 



P P 



where 



+ ^Ll{p^S{fp-'),r}C^nAp-'p-H{r - r{t))6{e - e{t)) 
-\p^A^SJ+{p-^J+{pp-^Crnrn5{r - r{t))5{6 - e{t)))} 



L\= de T + aw sin ^ + s cot 0, 



(26) 
(27) 



and S denotes -2 5"," (6*) for simplicity. 

For a source bounded in a finite range of r, it is convenient to rewrite Eq. (26) 
further as 



Timo. = M/ dte''^*-''"'^(*^A2 (A„„o + A^no + A^™o)^(r-r(t)) 

'J —OG 

+ {(^Tnni+A^miM(r-r(t))}^ + {A^,^,2<^(r-r(t))} J, (28) 



where 



-^nn 
-^rnn 

-^m m 
-^rnn 1 
-^m m 1 
-^m m 2 



27rA2 
2 



Cnnp-'p-'L+{p-''L+{p'S)}, 

iK 



Cm nP 



{LtS)(}^+p + ^) 



—a sin 9S—{p — p) 



-^=r-p ^Cjrin[LtS + iasiue{-p- p)S], 
^ zp-^pCrnmS [i^ + , 



2^' 

1 



P pCmmS- 



(29) 
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Inserting Eq. (28) into Eq. (20), we obtain Znjnui as 

^ ■/ (ite*'"*-'™'^(*)iy^„^, (30) 

J — oo 



where 



J pin j2 pin 

l^mn 1 T -^mm 1 J T ,0 ^mrn 2 • V'^-'-^ 



dr av^ 



r=r{t) 



In this paper, we focus on orbits which are either circular (with or witliout 
inchnation) or eccentric but confined on the equatorial plane. In cither case, the 
frequency spectrum of T£rnui becomes discrete. Accordingly, Zimu in Eq. (19) 
or (20) takes the form, 

Zimiv = ^ S{UJ - UJn)Zemui ■ (32) 
n 

Then, in particular, ^4 at r ^ oo is obtained from Eq. (8) as 

At infinity, V4 is related to the two independent modes of gravitational waves 
/i+ and /ix as 

i}i = \(h+-ih^). (34) 

Prom Eqs. (33) and (34), the luminosity averaged over t ;» At, where At is the 
characteristic time scale of the orbital motion (e.g., a period between the two 
consecutive apastrons), is given by 



/ dE\ _ ^ _ ^ ^dE\ 



\ dt /-^ AlTU^ /-^ V dt J tmn 

In the same way, the time-averaged angular momentum flux is given by 
I ^ 



dt I „^ AnOJ^ \ dt ) emn „^ UJn^ dt ) tmn' 

«,m,n t,m^n £^m,n 



2.2 Chandrasekhar-Sasaki-Nakamura transformation 

As seen from the asymptotic behaviors of the radial functions given in Eqs. (19) 
and (20), the Teukolsky equation is not in the form of a canonical wave equation 
near the horizon and infinity. Therefore, it is desirable to find a transformation 
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that brings the radial Teukolsky equation into the form of a standard wave 
equation. 

In the Schwarzschild case, Chandrasekhar found that the Teukolsky equation 

can be transformed to the Rcggc-Whcclcr equation, which has the standard 
form of a wave equation with solutions having regular asymptotic behaviors at 
horizon and infinity [19]. The Regge- Wheeler equation was originally derived as 
an equation governing the odd parity metric perturbation [61]. The existence 
of this transformation implies that the Regge- Wheeler equation can describe 
the even parity metric perturbation simultaneously, though the explicit relation 
of the Regge- Wheeler function obtained by the Chandrasekhar transformation 
with the actual metric perturbation variables has not been given in the literature 
yet. 

Later, Sasaki and Nakamura succeeded in generalizing the; Chandrasekhar 
transformation to the Kerr case [66, 67] . The Chandrasekhar-Sasaki-Nakamura 
transformation was originally introduced to make the potential in the radial 
equation short-ranged and to make the source term well-behaved at horizon 
and infinity. Since we are interested only in bound orbits, it is not necessary to 
perform this transformation. Nevertheless, because its flat space limit reduces to 
the standard radial wave equation in the Minkowski spacetimc, it is convenient 
to apply the transformation when dealing with the post-Minkowski or post- 
Newtonian expansion, at least at low orders of expansion. 

We transform the homogeneous Teukolsky equation to the Sasaki-Nakamura 
equation [66, 67], which is given by 

F{r)-^-U{r)\Xemu, = 0. (37) 
The function F(r) is given by 

= (38) 

V = co + ci/r + + cs/r^ + a/r'^, (39) 



where 
with 



Co = -UiojM + X{X + 2) - 12aijj{aLV - m), 

Cl = 8ia[3au} — X{au} — m)], 

C2 = -2AiaM{a(jO - m) + 12a^[l - 2{auj - mf], 

C3 = 24:ia^{au> - m) - 24Ma^, 

C4 = 12a''. (40) 



The function U{r) is given by 



U{r) = +G' + 4% - (41) 



13 



where 



G 

a 



2(r - M) 



A2 



rA 
(r2 + a2)2 

/3, 



2A(~iK 



2A\ 

- M . 

r / 



The relation between Rimui and Ximui is 



(42) 



(43) 



where Ximuj = Ximuj^/ (r^ + a'^Y^^. Conversely, we can express Ximu in terms 
of Remoj as 

X,„. = (r2 + a2)V'r-2j_J_ 



(44) 



where J_ = (d/rfr) - i{K/A). 

If wc set a = 0, this transformation reduces to the Chandrasekhar transfor- 
mation for the Schwarzschild black hole [19]. The explicit form of the transfor- 
mation is 

( d . ( d \ Rimuj 

= — — — \-tco]- 

A\dr* J A\dr* J 



(46) 



where co, defined in Eq. (40), reduces to cq = {£ - !)£{£ + !){£ + 2) - 12iMuj. 
In this case, the Sasaki-Nakamura equation (37) reduces to the Regge- Wheeler 
equation [61] which is given by 



d^ 



dr*"^ 



+ uj^- V{r) 



XeAr) = 0, 



where 



V{r) = 1 - 



2M\ fe{e + i) 6M 



(47) 



(48) 



As clear from the above form of the equation, the lowest order solutions are 

given by the spherical Bcsscl functions. Hence it is intuitively straightforward to 
apply the post-Newtonian expansion to it. Some useful techniques for the post- 
Newtonian expansion were developed for the Schwarzschild case by Poisson [56] 
and Sasaki [65]. 

The asymptotic behavior of the ingoing- wave solution X™ which corresponds 
to Eq. (14) is 



X, 



/Aiei piuir' i /linc p-iwr* 



rans^— ifcr* 



for r* 
for r* 



oo 



-oo. 



(49) 
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The coefficients A'"^ A'^^ and A*"*"^ are related to B'"", B'^^ and 5*^'*"% defined 
in Eq. (14), by 



r>inc 

^emw — 

r>ref 

^imui — 


/line 
4 ref 

Co 


(50) 
(51) 


Dtrans 


4 trans 


(52) 





where 

+{12iam - 16M + IGamMuj + 24iM^uj)r+ 
-4a^m^ - 12iamM + 8M^]. 

In the following sections, we present a method of post-Nowtonian expan- 
sion based on the above formalism in the case of the Schwarzschild background. 
In the Kerr case, although a post-Newtonian expansion method developed in 
previous work [69, 74] was based on the Sasaki-Nakamura equation, we will 
not present it in this paper. Instead, we present a different formalism, namely 
the one developed by Mano, Suzuki and Takasugi which allows us to solve the 
Teukolsky equation in a more systematic manner, albeit very mathematical [43]. 
The reason is that the equations in the Kerr case are already complicated enough 
even if one uses the Sasaki-Nakamura equation, so that there is not much ad- 
vantage in using it. In contrast, in the Schwarzschild case, it is much easier to 
obtain physical insight into the role of relativistic corrections if we deal with the 
Regge- Wheeler equation. 

3 Post-Newtonian expansion of the Regge- Wheeler 
equation 

In this section, we review a post- Newtonian expansion method for the Schwarzschild 
background, based on the Regge- Wheeler equation. We focus on the gravita- 
tional waves emitted to infinity, but not on those absorbed by the black hole. 
The black hole absorption is deferred to Section 4, in which we review the 
Mano-Suzuki- Takasugi method for solving the Teukolsky equation. 

Since we are interested in the waves emitted to infinity, as seen from Eq. (20), 
what we need is a method to evaluate the ingoing-wave Teukolsky function 
^Tmun O'' its counterpart of the Regge- Wheeler equation, X^^^, which are re- 
lated by Eq. (43). In addition, we assume w > whenever it is necessary 
throughout this section. Formulas and equations for a; < are obtained from 
the symmetry X]^^ = X'IL^_^. 
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3.1 Basic assumptions 

We consider the case of a test particle with mass /i in an orbit which is nearly 
circular around a black hole with mass M » /x. For a nearly circular orbit, 
say at r ~ ro, what we need to know is the behavior of B}(!^^ at r ~ tq. In 
addition, the contribution of uj to R™„^^ comes mainly from w ~ mO.^, where 
fi^ ~ (M/rg)^/^ is the orbital angular frequency. 

Thus, if we express the Regge- Wheeler equation (47) in terms of a non- 
dimensional variable z = lot, with a non-dimensional parameter e = 2Ma;, 
we are interested in the behavior of X^i^^[z) at z ~ wro ~ m{M/ro)^^^ ~ v 
with e ~ 2m(M/ro)^^^ ~ v^, where v = {M/roY^^ is the characteristic orbital 
velocity. The post-Newtonian expansion assumes that v is much smaller than 
the velocity of light; v <^ 1. Consequently, we have e <C u <C 1 under the 
post-Newtonian expansion. 

To obtain X™^^ (which we denote below by Xf for simplicity) under these 
assumptions, we find it convenient to rewrite the Regge- Wheeler equation in an 
alternative form. It is 



jf_ ld_ , ^(^+1)^ 
dz"^ z dz 



+ I' - -^)\^' = Tz [v^Tz ' ^'^'^^ 

where is a function related to X^ as 

X| = ze-'^'"(^-^)e^. (54) 

The ingoing-wave boundary condition of is derived from Eqs. (49) and (54) 
as 

/ ylf ^^e'^'^'z-^e"'^ + fe-'^'"^z-ie'("^+2^^"^^ for r* oo, 



1 d 



(53) 



(55) 



The above form of the Regge- Wheeler equation is used in the following subsec- 
tions. 

It should be noted that if we recover the gravitational constant G, we 
have e = 2GMu. Thus, the expansion in terms of e corresponds to the post- 
Minkowski expansion, and expanding the Regge- Wheeler equation with the as- 
sumption e <C 1 gives a set of iterative wave equations on the flat spacetime 
background. One of the most significant differences between the black hole per- 
turbation theory and any theory based on the flat spacetime background is the 
presence of the black hole horizon in the former case. Thus, if we naively expand 
the Regge- Wheeler equation with respect to e, the horizon boundary condition 
becomes unclear, since there is no horizon on the flat spacetime. To establish 
the boundary condition at the horizon, we need to treat the Regge- Wheeler 
equation near the horizon separately. We thus have to find a solution near the 
horizon and the solution obtained by the post-Minkowski expansion must be 
matched with it in the region where both solutions are valid. 

It may be of interest to note the difference between the matching used in the 
BDI approach for the post-Newtonian expansion [12, 15] and the matching used 
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here. In the BDI approach, the matching is done between the post-Minkowskian 
metric and the near-zone post-Newtonian metric. In our case, the matching is 
done between the post-Minkowskian gravitational field and the gravitational 
field near the black hole horizon. 



3.2 Horizon solution; z <^ 1 

In this subsection, we first consider the solution near the horizon, which we call 
the horizon solution, based on [59]. To do so, we assume 2: <C 1 and treat e as 
a small number, but leave the ratio z/e arbitrary. We change the independent 
variable to a; = 1 — z/e and the wave function to 



Itr; 

Note that the horizon corresponds to x = 0. We then have 



x{x - l)Z" + [2(3 - ie)x - (1 - 2ie)] Z' 

+ [6 - t{£ -hi) - 5ie + e^{2-2,x + x^)] Z = 0, (57) 

where a prime denotes differentiation with respect to x. We look for a solution 
which is regular at x = 0. 

First, we consider the lowest order solution by setting e = in Eq. (57). 
The boundary condition (55) requires that Z = 1 at x = 0. The solution which 
satisfies the boundary condition is 



e-2 



n=0 



Thus, the lowest order solution is a polynomial of order i? — 2ina; = l — z/e. 

Next, we consider the solution accurate to 0(e). We neglect the terms of 
0{e^) in Eq. (57). Then its takes the form of a hypergeometric equation, 

x{x - l)Z" +[{a + h+ l)x -c]Z' + abZ = 0, (59) 

with parameters 

a = -{i -2) -ie + Oie"^), 

b = £ + 3 - ie + Oie"^), (60) 
c = 1 - 2ie. 

The two linearly independent solutions are F{a, b; c; x) and x^^'-F{a + 1 — c,b + 
1 — c;2 — c;x), where F is the hypergeometric function. However, only the first 
solution is regular at a; = 0. We therefore obtain 

^i{z)=Af^'''e-^e''^^'''-^^ (^)'F(a,6;c;l- ^) . (61) 
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The above solution must be matched with the solution obtained from the 
post-Minkowski expansion of Eq. (53), which we call the outer solution, in a 
region both solutions are valid. It is the region where the post-Newtonian 
expansion is applied, i.e., the region e ^ z <C 1. For this purpose, we rewrite 
Eq. (61) as (see, e.g., Eq. (15.3.8) of [7]) 



^ trans —l_ie(ln e — 1) 



+ 



'^+i<^ r(c)r(6 - a) 



ej r(6)r(c 

1+^^ r(c)r(a - h) 



r(a)r(c-6) 



a) 

F{b, c — a;b 



F{a, c — b;a- 



6 + 1;-) 

z 



1;-) 

z 



(62) 



This naturally allows the expansion in e/z. It should be noted that the second 
term in the square brackets of the above expression is meaningless as it is, since 
the factor r(a — b) diverges for integer i. So, when evaluating the second term, 
we first have to extend £ to a non-integer number. Then, only after expanding 
it in terms of e, we should take the limit of an integer £. One then finds that this 
procedure gives rise to an additional factor of 0(e). For e/z <C 1, it therefore 
becomes ©(e^^"*"^) higher in e than the first term. Then we obtain 



-,i r 



2)!(£ + 2)!e'^+i 
(£-2)(£ + 2)e 
2i z 



1 -|- ie{ai + In 2;) 
+ 0(e2) , 



where 

a^ = 27 + V^(^-l)+V(^ + 3)-l, 
and ^(n) is the digamma function, 

n-l 

V;(n) = -7 + ^ 
fe=i 



(63) 
(64) 

(65) 



and 7 ~ 0.57721 is the Euler constant. 

As we will see below, the above solution is accurate enough to determine the 
boundary condition of the outer solution up to 6PN order of expansion. 



3.3 Outer solution; e <C 1 

We now solve Eq. (53) in the limit e <C 1, i.e., by applying the post-Minkowski 
expansion to it. In this subsection, we consider the solution to 0(e). Then we 
match the solution to the horizon solution given by Eq. (63) at e <C ^ <IC 1. 
By setting 

00 

n=0 

each ^^"^ {z) is found to satisfy 



(66) 



£_ 2d_ 
dz'^ z dz 



1 



£{£ + !) 



in) 



dz 



]_d_ 
z^ dz 



(67) 
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Equation (67) is an inhomogcncous spherical Bcsscl equation. It is the simphc- 
ity of this equation which motivated the introduction of the auxihary function 
[65]. 

The zeroth-order solution, satisfies the homogeneous spherical Bessel 
equation, and must be a linear combination of the spherical Bessel functions of 
the first and second kinds, ji{z) and n^(z). Here, we demand the compatibility 



with the horizon solution (63). Since jt{,z) ~ and rn{z) 



-e-i 



ni{z) does 



not match with the horizon solution at the leading order of e. We therefore have 

£,f\z) = afM{z). (68) 

The constant a^P represents the overall normalization of the solution. Since it 

can be chosen arbitrary, wo sot = 1 below. 

The procedure to obtain ^^^■'(z) was described in detail in [65]. Using the 
Green function G{z, z') = je{z<)ni{zy), Eq. (67) may be put into an indefinite 
integral form, 



f 



dzz'^e '"je 



f 



dzz e "rii 



lie-z'4-'\z)y 



(69) 



The calculation is tedious but straightforward. All the necessary formulas to 
obtain ^^"^ for n < 2 are given in Appendix of [65] or Appendix D of [46] . Using 
those formulas, we have for n = 1, 



6'^ = a^'ji + P^^'n, 



(1), 



t-2 



( e-A 



21- 1 



V 2^(2^+1) 



Je-i 



/I 1 \ 
+Refijo + — + — -r Re,mjm - ^D"^^ + ije In; 



(70) 



Here, and Rt m are functions defined as follows. The function D^-' is given 
by 

(71) 



= \ be Si{2z) - m {C\{2z) - 7 - In 2^)] , 



where C\{x) = — dtcost/t and Si(a;) = dtsint/t. The function Rm,k is 
defined by Rm,k = z'^{nmjk — jmnk)- It is a polynomial in inverse powers of z 
given by 



Rm,k — 

for m> k, and 



r=0 



! r(rn - A; - 2r)r (A; + | + r) 



m—k—l—2r 



(72) 
(73) 
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for m < k. 

Here, we again perform the matching with the horizon solution (63). It 
should be noted that given by Eq. (70) is regular in the limit z ^ Q except 
for the term ff^^nt. By examining the asymptotic behavior of Eq. (70) at 
2 <C 1, we find /3^^^ = 0, i.e., the solution is regular at 2; = 0. As for Oi'p , it 

only contributes to the renormalization of Hence we set a^^^ = 0. Thus, 
the transmission amplitude Af'^'^^ is determined to 0(e) as 

It may be noted that this explicit expression for Af'^''^ is unnecessary for the 
evaluation of gravitational waves at infinity. It is relevant only for the evaluation 
of the black hole absorption. 



3.4 More on the inner boundciry condition of the outer 
solution 

In this subsection, wc discuss the inner boundary condition of the outer solution 
in more detail. As we have seen in the previous subsection, the boundary 
condition of is that it is regular a.t z ^ 0, at least to 0(e). While in the full 
non-linear level, the horizon boundary is at z e. We therefore investigate to 
what order in e the condition of regularity at 2: = can be applied. 

Let us consider the general form of the horizon solution. With x = 1 — z/e, 
it is expanded in the form, 

= > (x) + e ei'^ (x) + {x) + ---. (75) 

The lowest order solution ^j°^(a;) is given by the polynomial (58). Apart from 
the common overall factor, it is schematically expressed as 



l + ci - -i h ce-2 ( - 

z \z 



(76) 



Thus, ^j^-*^ does not have a term matched with ng, but it matches with j^. 
We have = z^e~^ ~ e~^ji. A term that matches with first appears 
from ■ This can be seen from the horizon solution valid to 0(e), Eq. (62). 
The second term in the square brackets of it produces a term e{z/e)~^~^ = 
^e+2^-e-i ^ c^+'^ni. This term therefore becomes 0(e^^+^/2;^^+^) higher than 
the lowest order term Since £ > 2, this effect first appears at 0(e®) 

in the post-Minkowski expansion, while it first appears at 0{v^^) in the post- 
Newtonian expansion if we note that e = 0{v^) and z = 0{v). This implies, 
in particular, that if we are interested in the gravitational waves emitted to 
infinity, we may simply impose the regularity at z = as the inner boundary 
condition of the outer solution for the calculation up to 6PN order beyond the 
quadrupole formula. 
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The above fact that a non-trivial boundary condition due to the presence 
of the black hole horizon appears at ©(e^^"''-^) in the post-Minkowski expansion 
can be more easily seen as follows. Since je = 0{z^) as z ^ 0, we have Xd — > 
0(e^+i)e-^^*, or A*™'' = ©(e'^+i), where z* = z + eln(z - e). On the other 
hand, from the asymptotic behavior of at 2; = oo, the coefficients A™^ and 
A^^^ must be of order unity. Then using the Wronskian argument, we find 

j 4 trans 1 2 

Thus, we immediately see that a non-trivial boundary condition appears at 
0(e2^+2). 

It is also useful to keep in mind the above fact when we solve for under the 

post-Minkowski expansion. It implies that we may choose a phase so that A^^^ 
and are complex conjugate to each other, to 0(e^^+^). With this choice, 
the imaginary part of Xg, which reflects the boundary condition at horizon, do 
not appear until 0(e^^+^) because the Regge- Wheeler equation is real. Then 

recalling the relation of to X^, Eq. (54), Im ^C^"^^ for a given n < 2^ + 1 is 

completely determined in terms of Re (it^^^ for r < n — 1. That is, we may 
focus on solving only the real part of Eq. (67). 

3.5 Structure of the ingoing-wave function to O(e^) 

With the boundary condition discussed in the previous section, we can integrate 

the ingoing- wave Rcggc- Wheeler function iteratively to higher orders of e in the 
post-Minkowskian expansion, e <C 1. This was carried out in [65] to O(e^) and 
in [77] to O(e^) (See [46] for details). Here, we do not recapitulate the details 
of the calculation since it is already quite involved at O(e^), with much less 
space for physical intuition. Instead, we describe the general properties of the 
ingoing- wave function to O(e^). 

As discussed in the previous section, the ingoing wave Regge- Wheeler func- 
tion Xi can be made real up to 0(e^^+^), or to O(e^) of the post-Minkowski 
expansion if we recall i > 2. Choosing the phase of Xe in this way, let us ex- 
plicitly write down the expressions of Im(^^"-') (n = 1.2) in terms of Re(^^'"^) 
(to < n — 1). We decompose the real and imaginary parts of as 

^W=/W+i5W. (78) 

Inserting this expression into Eq. (54) and expanding the result with respect to 
e, and noting f^^^ — ji and gf^ = 0, we find 

X, = e--'"(^-^)z(j, + e(/(^) + z^f ) + e^fP + igf^) + ■ ■ ■) 
= ^ [k + e/i'^ + (/i^) + 5« In. - ij,(ln^)^) + • 

+iz f e(<?f ) - M In z) + ( gf ^ + -jt - \nz\+--\. (79) 
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Hence, we have 

gf^=je\nz = j> + /f hn^, • • • . (80) 

We thus have the post-Minkowski expansion of as 

oo 
n=0 

Xf^=zje, xW=./«, =.(/f +ii,(ln.)^),---.(81) 

Now, let us consider the asymptotic behavior of at <C 1. As we know 
that and ^ arc regular at z = 0, it is readily obtained by simply assuming 
Taylor expansion forms for them (including possible In 2; terms), inserting them 
to Eq. (67), and comparing the terms of the same order on both sides of the 
equation. Wc denote the right-hand side of Eq. (67) by Sj"^ . 

For n = 1, we have 

4 + z2 ■ 



r2 



-Je. 



0{z) for ^ = 2, 
0{z^-^) for ^ > 3. 



(82) 



Inserting this to Eq. (67) with n = 1 we find 



Rp _ _ i ^(z^) for ^ = 2, , . 

Re {^t }-ft -\ o{z'-') for £ > 3. 

Of course, this behavior is consistent with the full post-Minkowski solution given 
in Eq. (70). 

For n = 2, we then have 



^ , V ^ . , , J 0(z^-2) for f = 2,3, 
--(,,ln.) --^,,lnz+|^J^,_,j ^^^^^^1 (84) 



This gives 

R,/.(2A _ .(2)_ rO(z^) + 0(z^)lnz-ij,(lnz)2 for £ = 2, 3, 



Note that the Inz terms in (84) arising from g^''''^ gives the (Inz)^ term in /j^'' 

that just cancels the j^(ln 2)^/2 term of xj^^ in Eq. (81). 

Inserting Eqs. (83) and (85) into the relevant expressions in Eq. (81), we 
find 

X2 = z^[0{l) + eO{z) + e''{0{l)+0{l)lnz} + ■■■], 
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Xa = [Oiz] + eOil) + e^O{z) + Oiz) In z} + ■■■], (86) 
Xe = z^[0{z^-^) + €0{z^-^) + e''{0{z^-^) + 0{z^-^)lnz} + ---] {£>4). 

Note that, for £ = 2 and 3, the leading behavior of X^"' at n = ^ — 1 is more 
regular than the naively expected behavior, ^ which propagates to the 

consecutive higher order terms in e. This behavior seems to hold for general i, 
but we do not know a physical explanation of it. 

Given a post-Newtonian order to which we want to calculate, by setting 
z = 0{v) and e = O(u^), the above asymptotic behaviors tell us the highest 
order of X^"^ we need. We also see the presence of In^; terms in xf'\ The 
logarithmic terms appear as a consequence of the mathematical structure of the 
Regge- Wheeler equation at z <C 1. The simple power series expansion of X^"^ in 
terms of z breaks down at O(e^), and we have to add logarithmic terms to obtain 
the solution. These logarithmic terms will give rise to Inu terms in the wave- 
form and luminosity formulas at infinity, beginning at 0{v^) [72, 73]. It is not 
easy to explain physically how these In v terms appear. But the above analysis 
suggests that the In v terms in the luminosity originate from some spatially local 
curvature effects in the near zone. 

Now we turn to the asymptotic behavior at = oo. For this purpose, let 
the asymptotic form of f\^^ be 

./I"' ^ Pt'^Ji + Qf^T^i as z -> 00 . (87) 

Noting Eq. (81) and the equality e"*'^"^*^^"*^) = e~*^*e*^, the asymptotic form of 
Xi is expressed as 



+^ { {PP - 0« \nz}+i (Qf + ) In .) } 



AT = -i^+'^e-'^^''^ 



Note that 



r — 2M' 

ojr* = w{ r + 2Mln— — — \ = z* - elne, 
' 2M ' 



because of our definition of z*, z* = z + e + \n{z — e). The phase factor e~'^'"^ 
of AT originates from this definition, but it represents a physical phase shift 
due to wave propagation on the curved background. 

As one may immediately notice, the above expression for AT contains Inz- 
dependent terms. Since AT should be constant, P^^^ and Q^"^ should contain 
appropriate In z-dependent terms which exactly cancel the In z-dependent terms 
in the formula (88). To be explicit, we must have 

p(i) _ „W 

= 5«-lnz, 
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Q 



5(2) 

e 

(2) 



J2) 



(89) 



where p^' and gj"^ are constants. These relations can be used to check the 
consistency of the solution obtained by integration. In terms of p^"' and 



An) 



A™'^ is expressed as 



(90) 



Note that the above form of A^^^ imphes that the so-called tail of radiation, 
which is due to the curvature scattering of waves, will contain In w terms as phase 
shifts in the waveform, but will not give rise to such terms in the luminosity 
formula. This supports our previous argument on the origin of the Inv terms 
in the luminosity. That is, it is not due to the wave propagation effect but due 
to some near-zone curvature effect. 



4 Analytic solutions of the homogeneous Teukol- 
sky equation by means of the series expansion 
of special functions 

In this section, we review a method developed by Mano, Suzuki and Taka- 
sugi [43] who found analytic expressions of the solutions of the homogeneous 
Teukolsky equation. In this method, the exact solutions of the radial Teukolsky 
equation (9) are expressed in two kinds of series expansions. One is given by 
a series of hypergeometric functions and the other by a series of the Coulomb 
wave functions. The former is convergent at horizon and the latter at infinity. 
The matching of these two sohitions arc done exactly in the overlapping region 
of convergence. They also found that the series expansions are naturally related 
to the low frequency expansion. Properties of the analytic solutions were stud- 
ied in detail in [45]. Thus, the formalism is quite powerful when dealing with 
the post-Newtonian expansion, especially at higher orders. 

In many cases when we study the perturbation of a Kerr black hole, it is 
more convenient to use the Sasaki-Nakamura equation since it has the form of 
a standard wave equation, similar to the Regge- Wheeler equation. However, it 
is not quite suited for investigating analytic properties of the solution near the 
horizon. In contrast, the Mano- Suzuki- Takasugi (MST) formalism allows us to 
investigate analytic properties of the solution near the horizon systematically. 
Hence, it can be used to compute the higher order post-Newtonian terms of the 
gravitational waves absorbed into a rotating black hole. 

We also note that this method is the only existing method that can be used 
to calculate the gravitational waves emitted to infinity to an arbitrarily high 
post-Newtonian order in principle. 
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4.1 Angulcir eigenvalue 



The solutions of the angular equation (10) that reduce to the spin-weighted 
spherical harmonics in the limit aw ^ are called the spin-weighted spheroidal 
harmonics. They are the eigen functions of Eq. (10) with A being the eigenvalues. 
The eigenvalues A are necessary for discussions of the radial Teukolsky equation. 
For general spin weight s, the spin weighted spheroidal harmonics obeys 

sin 6* ^6* I dQ) sin^ 6* 

— 2(Ui)s cos ^ + s + 2mau! + A sSe-m = 0. (91) 

In the post-Newtonian expansion, the parameter aw is assumed to be small. 
Then, it is straightforward to obtain a spheroidal harmonic sSem of spin-weight 
s and its eigenvalue A perturbatively by the standard method [60, 74, 69]. 

It is also possible to obtain the spheroidal harmonics by expansion in terms 
of the Jacobi functions [29]. In this method, if we calculate numerically, we can 
obtain them and their eigenvalues for arbitrary value of acu. 

Here we only show an analytic formula for the eigenvalue A accurate to 
O ((aw)^) , which is needed for the calculation of the radial functions. It is given 
by 

X = Xo+ acuXi + a^iJ'Xi + 0{{awf) (92) 



where 



Ao =^(£+l)-s(s + l), 
Ai = -2m ( 1 



with 



As = 1 + {n{i + 1) - n{i) - 1), (93) 



^ ^ (2f- 1)^3(2^+1) • y^'^) 



4.2 Horizon solution in series of hypergeometric functions 

As in Sec. 3, we focus on the ingoing-wave function of the radial Teukolsky 
equation (9). Since the analysis below is applicable to any spin, |s| = 0, 1/2, 1, 
3/2 and 2, we do not specify it except when it is needed. Also, the analysis is 
not restricted to the case aw <C 1 unless so stated explicitly. For general spin 
weight s, the homogeneous Teukolsky equation is given by 

As before, taking account of the symmetry Riimui = Re-m-u, we may assume 
e = 2Mu > if necessary. 
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The Tciikolsky equation has two regular singularities at r = r±, and one 
irregular singularity at r = oo. This implies that it cannot be represented in the 
form of a single hypergeometric equation. However, if we focus on the solution 
near the horizon, it may be approximated by a hypergeometric equation. This 
motivates us to consider the solution expressed in terms of a series of hyperge- 
ometric functions. 

We define the independent variable x in place of z (= wr) as 



en 

where 



(96) 



M 

For later convenience, we also introduce t — {e — mq)/K and e± = (e ± t)/2. 
Taking into account the strTicture of the singularities at r = r±, we put the 
ingoing- wave Teukolsky function R^^^ as 

P^Lu. = e"'^^(-.T)-^-*(^+^)/2(l - xf '-^^'^p,^{x). (97) 

Then the radial Teukolsky equation becomes 

a;(l — x)piJ' + [l — s — ie ~ IT — [2 — 2ir)a;]pin' 
+ [iT{l -iT) + \ + s{s + l)]pi„ 
= 2ieK[—x{\ — x)pin' + {1 — s + ie — iT)xpin] 

+[e'^ - i€K{l - 2s)]pi„ , (98) 

where a prime denotes d/dx. The left-hand-side of Eq. (98) is in the form of a 
hypergeometric equation. In the limit e ^ 0, noting Eq. (92), we find that a 
solution which is finite at a; = is given by 

Pi„(e ^ 0) = F{-1 -iT,£+l-iT,l- s- ir, x). (99) 

For a general value of e, Eq. (98) suggests that a solution may be expanded 
in a series of hypergeometric functions with e being a kind of an expansion 
parameter. This idea was extensively developed by Leaver [42]. Leaver ob- 
tained solutions of the Teukolsky equation expressed in a series of the Coulomb 
wave functions. The MST formalism is an elegant reformulation of the one by 
Leaver [42]. 

The essential point is to introduce the so-called renormalized angular mo- 
mentum u which is a generalization of £ to a non-integer value such that the 
Teukolsky equation admits a solution in a convergent series of hypergeometric 
functions. Namely, we add the term [i/(i/ -t- 1) — A — s{s + l)]pin to both sides 
of Eq. (98) to rewrite it as 

x{l — x)pin" + [I — s — ie — ir — {2 — 2iT)x]pin' 

+ [iT{l - it) + V{V + l)]pi„ 

= 2ieK[—x{l — x)pi^ -|- (1 — s -I- ie — ir)a:pin] 

+[v{u + 1) - A - s(s -h 1) + - ieK(l - 2s)]pin ■ (100) 



26 



Of course, no modification is done to the original equation, and v is just an 
irrelevant parameter at this stage. A trick is to consider the right-hand-side of 
the above equation as a perturbation, and look for a formal solution specified 
by the index v in a scries expansion form. Then, only after wc obtain the 
formal solution, we require that the series should converge and this requirement 
determines the value of v. Note that, if we take the limit e ^ 0, we must have 
V ^ i (ox V ^ —t — 1), to assure \y{y -|- 1) — A — s{s -|- 1)] ^ and to recover 
the solution (99) 

Let us denote the formal solution specified by a value of by pj^. We express 
it in the series form, 



(101) 



Pn+i^{x) = F{n + v + 1 — iT,—n — v — ir;! — s — ie — iT;x). (102) 
Here, the hypergeometric functions Pn+v{x) satisfy the recurrence relations [43], 

{n + v + l — s — ie){n + v + 1 — ir) 



XPn+v — 



+ 



1 + 



2(n + ;/-hl)(2r7, + 2;y + l) 
iT{s + ie) 



-Pn+v+l 



{n + v){n + v+l) 



Pn+v 



{n + v + s + ie){n + y + iT) 
2{n + v){2n + 2u +1) 



-Pn+u-l, 



(103) 



{n + V + ir) {n + V + 1 — ir) {n + u +1 



ze) 



+ i(s + ie) 



2{n + v +l){2n + 2v +1) 

zr(l — it) 



Pn+u+l 



(n + z/)(n + z/+ 1) 



{n + v + 1 — iT){n + + iT){n + + s + ie) 
2{n + u){2n + 2u + 1) 



p„+^_i,(104) 



Inserting the series (101) into Eq. (100) and using the above recurrence relations, 
we obtain a three-term recurrence relation among the expansion coefficients a„. 
It is given by 

a'^ttn+l + Pnttn + Jn^n-l = 0, (105) 

where 



a„ = 



0" 



In 



i€K{n -|-i'-|-l-|-s-|-ie)(n-|-i'-|-l-|-s — ie){n + 1/ + l + iT) 

{n + v + l){2n + 2iy + 3) 
-A - s(s +l) + {n + u){n + u + 1) + e'^ + e{e - mq) 

e(e — mq){s'^ + f^) 
(n-|-z/)(n + /y + l)' 

icK(ii + — -s- + /'()(" + ^ ^ '< )(" + '-' ^ 't) 
{n + i'){2n + 2u - 1) " 



(106) 
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The convergence of the series (101) is determined by the asymptotic behaviors 
of the coefficients aj^ at n ^ ±00. We thus discuss properties of the three- term 
recurrence relation (105), and the role of the parameter v in detail. 

The general solution of the recurrence relation (105) is expressed in terms 
of two linearly independent solutions {/jl^^} and {/n^^} (n = ±1, ±2, • • •). Ac- 
cording to the theory of three-term recurrence relations ([35], p. 31) when there 
exists a pair of solutions which satisfy 

/■(I) / .(1) \ 

hm ^ = hm ^ = , (107) 

then the solution {/rl^^} is called minimal as n — > 00 (n ^ — 00). Any non- 
minimal solution is called dominant. The minimal solution (either as n ^ 00 
or as n — > —00) is determined uniquely up to an overall normalization factor. 

The three-term recurrence relation is closely related to continued fractions. 
We introduce 

Rn^^, Ln^^. (108) 
dn-l dn+l 

We can express -R„ and L„ in terms of continued fractions as 



Rn — — 



7n <^nln+l ««+l7n+2 



(109) 



(110) 



These expressions for Rn and i„ are valid if the respective continued fractions 

converge. It is proved ([35], p. 31) that the continued fraction (109) converges if 
and only if the recurrence relation (105) possesses a minimal solution as n ^ 00, 
and the same for the continued fraction (110) as n — > —00. 

Analysis of the asymptotic behavior of (105) shows that, as long as v is 
finite, there exists a set of two independent solutions which behave as (e.g., [35], 
p.35), 

(1) (2) o- 

,. ail' ten ait' 2i 

hm n-rp- = — , lim — ^ = — , (111) 
and another set of two independent solutions which behave as 

lim n . = , lim — tttt— = . (112) 

Thus {an''} is minimal as n — > 00 and ^} is minimal as n ^ —00. 



28 



Since the recurrence relation (105) possesses minimal solutions as n ^ ±00, 
the continued fractions on the right-hand-sides of Eqs. (109) and (110) converge 
for a„ = ttn^ and a„ = bn\ In general, however, and bn^ do not coincide. 
Here, we use the freedom of v to obtain a consistent solution. Let {/^} be a 
sequence which is minimal for both n ^ ±00. We then have expressions for 
fn/fn-i ^^'-^ fn/fn+1 terms of Continued fractions as 

5 _ fn _ In (^^'^\ 

Jn-1 Pn Pn+l Pn+2 



fn OCn OLn-\ln Ot-n-^Hn 



fn+1 a^n- P'n-1- Pn-2- 



(114) 



This implies 



RnLn-l = l. (115) 

Thus, if we choose v such that it satisfies the implicit equation for v, Eq. (115), 
for a certain n, we obtain a unique minimal solution {fn} which is valid over 
the entire range of n, —00 < n < 00, that is 

linin^ = ^, hm n^ = -^. (116) 

Jn—1 Jn+1 

Note that if Eq. (115) for a certain value of n is satisfied, it is automatically 
satisfied for any other value of n. 

The minimal solution is also important for the convergence of the series 
(101). For the minimal solution {f^}, together with the properties of the hy- 
pergeometric functions Pn+,y for large |n|, we find 

i^) fn-lVn + v-lix) 

iim n — — = — hm n — — 

= !^[l-2x + ((l-2x)2-l)V2]. (117) 

Thus the series of hypergeomctric functions (101) converges for all x in the range 
> a; > —00 (in fact, for all complex values of x except at \x\ = 00), provided 
the coefficients are given by the minimal solution. 

Instead of Eq. (115), wc may consider an equivalent but practically more con- 
venient form of an equation that determines the value of v. Dividing Eq. (105) 
by a„, we find 

/?^ + <ii„+i+7;:j:n-i=0, (118) 

where Rn+i and are those given by the continued fractions (113) and 

(114), respectively. Although the value of n in this equation is arbitrary, it is 
convenient to set n = to solve for y. 

For later use, we need a series expression for R™ with better convergence 
properties at large |a;|. Using analytic properties of hypergeometric functions, 
we have 
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where 



■0 



e' 



OO 



r(l - s - ie - iT)r{2n + 2i/ + 1) 



X 



n 



■v 



r(n + z/ + 1 - iT)r(n + 1/ + 1 - s - ie) 



n=— oo 



X (1 — a;)"'F(— n — — ir, — n — — s — ie; — 2n — 2z/; 



1 



).(120) 



l-x 



This expression expUcitly exhibits the symmetry of -Rin under the interchange 
of u and —v — 1. This is a result of the fact that u{u + 1) is invariant under 
the interchange i/ ^ — 1. Accordingly, the recurrence relation (105) has the 
structure that fZn~^ satisfies the same recurrence relation as 

Finally, we note that if is a solution of Eq. (115) or (118), v + k with an 
arbitrary integer k is also a solution, since v appears only in the combination of 
n-\-v. Thus, Eq. (115) or (118) contains infinite number of roots. However, not 
all of these can be used to express a solution we want. As noted in the earlier 
part of this subsection, in order to reproduce the solution in the limit e 0, 
Eq. (99), we must have v ^ i {or v ^ —i — 1 hy symmetry). We thus impose 
a constraint on v such that it must continuously approach ^ as e — > 0. 

4.3 Outer solution in series of Coulomb wave functions 

The solution in series of hypergeometric functions discussed in the previous 
subsection is convergent at any finite value of r. However, it does not converge 
at infinity, and hence the asymptotic amplitudes, B™'^ and B^^^ , cannot be 
determined from it. To determine the asymptotic amplitudes, it is necessary to 
construct a solution which is valid at infinity and match the two solutions in a 
region where both solutions converge. The solution convergent at infinity was 
obtained by Leaver in series of Coulomb wave functions [42]. In this subsection, 
we review Leaver's solution based on [45] . 

In this subsection again, by noting the symmetry Rimu} = Re-m-u}, we 
assume w > without loss of generality. 

First, we define a variable z = u){r — r_) = eK{l — x). Let us denote a 
Teukolsky function by Rc- We introduce a function f{z) by 




(121) 



Then the Teukolsky equation becomes 



z^f" + [z^ + (2e + 2is)z - A - s(s + 1)]/ 
= eKz{f" + /) + eK(s - 1 + 2ie)f 

— ■^■(fi; — i{e — mq)){s — 1 + ie)/ 
z 

+[-2e^ + emq + K(e^ + ies)]f . 



(122) 
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Wc sec that the right-hand-side is cxpUcitly of 0(e) and the left-hand-side is 
in the form of the Coulomb wave equation. Therefore, in the limit e ^ 0, we 
obtain a solution, 

f{z) = Fi{-is-e,z), (123) 
where -FL(r?, z) is a Coulomb wave function given by 

Fl{v, z) = e--2^5^+i^^^Li_!^$(L + 1 _ 2i + 2; 2iz), (124) 

and $ is the regular confluent hypergeometric function function (see [7] , section 
13) which is regular at z = 0. 

In the same; spirit as in the previous subsection, wc introduce the renormal- 
ized angular momentum v. That is, we add [A -|- s(s + 1) — i'{i'+l)]f{z) to both 
sides of Eq. (122) to rewrite it as 

z^f" + [z^ + (2e + 2is)z - v{v + 1)]/ 

= eKz{f" + f) + eK{s - 1 + 2ie+)/' - 4(k - i{e - mq)){s - 1 + ie)/ 

+ [-v{v + 1) + A -I- s(s + 1) - 2e2 + emq + K{e^ + i€s)]f . (125) 

We denote the formal solution specified by the index v by fu{z), and expand it 
in terms of the Coulomb wave functions as 

n= — oo ^ ' 

where (a)„ = r(a + n)/r(a). Then, using the recurrence relations among 

1 _ (n + + 1 + s - k) ^ is^e 

^ ~ (n + z/+l)(2n + 2i.+ l) "+"+^^ (n + !/)(n + z/+l) 

(n + v-s + ie) 

+ (n + ^)(2n + 2i. + l)^"+''-^' ^^^^^ 
, _ (n-Ft/)(n + iy + l + s-ze) ^ is + e ^ 

^"+- - + J, + i)(2n + 2z/ + 1) "+'^+' + („ + z,)(n + + 1) '^"+'^ 

(n + :^-H)(n-h;/-s-Hie) ^ 
(n + i')[2n -I- 2i^ -f 1) 

we can derive the recurrence relation among 6„. The result turns out to be 
identical to the one given by Eq. (105) for a„. We mention that the extra factor 
(z^+1 + .s — «e)„/(i^+l — s + «e)„ in Eq. (126) is introduced to make the recurrence 
relation exactly identical to Eq. (105). 

The fact that we have the same recurrence relation as Eq. (105) implies that 
if wc choose the parameter v in Eq. (126) to be the same as the one given by 
a solution of Eq. (115) or (118), the sequence {/^} is also the solution for {6„} 
which is minimal for both n — > ±c». Let us set 

^""^ ^ (i. + l-s + ie)„-^"- 
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By choosing v as stated above, we have the asymptotic value for the ratio of 
two successive terms of as 

lim n = hm n = — . (129) 



Using an asymptotic property of the Coulomb wave functions, we have 

We thus find that the series (126) converges at z > eK or equivalently r > r-|_. 

The fact that we can use the same i/ as in the case of hypcrgcometric func- 
tions to obtain the convergence of the series of the Coulomb wave functions is 
crucial to match the horizon and outer solutions. 

Here, we note an analytic property of the confluent hypergeometric function 
([28],p.259), 

$(a, c; a;) = e""^^(a, c; x) + IMe-(»--)e^^'(c - a, c; -x), (131) 

1 (C — ffl) 1 (o) 

where ^ is the irregular confluent hypergeometric function, and Im(a;) > is 
assumed. Using this with the identifications, 

a = n + iy + l — s + ie, 
c = 2(n + z^ + l), 
X = 2iz , 

we can rewrite R'^ (for w > 0) as 

i?^ = i?^+i?^, (132) 



where 



T{u + l + s-ie) ^ ' 



X 

n= — oo 



J2 «"/n(2-S)"*(" + 1/ + 1 - s + ie, 2n + 2i/ + 2; 2iz), 

(133) 



x^{n + i^ + l + s-ie,2n + 2iy + 2; -2iz). (134) 
By noting an asymptotic behavior of \E'(a, c;x) at large \x\, 

^{a,c;x) ^ x'"" as |a;| ^ oo, (135) 
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we find 



R+ = A:;^-ie-'(^+^'"^\ (136) 
R". = A!:^-i-2«e'(^+''"^\ (137) 



where 

+00 

T{u + l + s- ie) 



+ V lJ-\-^-i^s — iA ^ ■'^ ^ ' 



n— — 00 

+00 



A" = 2-i-«+-e-5^(^+i+^)e-^^ V + ^ ''j" r.(139) 

n— — 00 

We can see that the functions R'^ and R'^ are incoming- wave and outgoing- wave 
solutions at infinity, respectively. In particular, we have the upgoing solution, 
defined for s = —2 by the asymptotic behavior (15), expressed in terms of a 
series of Coulomb wave functions as 

i?"P = R'i . (140) 

4.4 Matching of horizon and outer solutions 

Now, we match the two types of solutions Rq and Rq. Note that both of them 
are convergent in a very large region of r, namely for en < z < 00. We sec that 
both solutions behave a z'^ multiplied by a single-valued function of z for large 
\z\. Thus, the analytic properties of Rq and Rq are the same, which implies 
that these two are identical up to a constant multiple. Therefore, we set 

= K^Rh- (141) 

In the region ck < 2 < 00, we may expand both solutions in powers of z 
except for analytically non-trivial factors. We have 

^ \ —s—i€+ 00 00 

Z \ « ^ , ^ 



= e""'e-'^{eK)-'-"'+z''+'''+ [ ij Z^Z]*^".. 

n= — 00 j=0 



-s — iej^ 00 00 



neK)-''-^^+r+-+ (--1) Yl 

k= — oo n=k 

(142) 

i€-\- 00 00 



^ n= — ooj=0 

(\ —s—ie-\- 00 k 
— -1] J2 J2 ^ 



fe=— 00 n= — 00 

(143) 
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where 



r(l - s - 2ie+)r(2n + 2v+l) 
T{n + V + ir)r(n + v + s- 



{-n-v- it) J {-n-v-s- ie)j ^^^^^-n+, 



ze) 



r{n + v + I - s + ie) ju + I + s - ie)n 
T{2n + 2iy + 2) {u + I - s + ie)n 



{n + u+l-s + ie)j 

~ Jn 



(144) 



(145) 



{2n + 2u + 2)j{j\) 

Then, by comparing each integer power of z in the summation, in the region 
CK <^ z < oo, and using the formula r(^)r(l — z) = n/ sinTr^, we find 

-1 



^(eK)"-"2- 
e*''=(2eK) 




^2-*i'^r(l -s- 2ie+)T{r + 2iy + 2) 



r{r + u + 1 - s + ie)r{r + z/ + 1 + zr)r(r + + I + s + ie) 




l)r(n 



1 + s + ie) r(n + u + l + iT) 



T{n + V + 1 - s - ie)T{n + u + 1 - ir) 
-1 



(i/ + 1 + s - ie)n ,1, 

J n 



(r - n)!(r + 2i/ + 2)„ (i/ + 1 - s + ie) 



fnj 

(146) 



where r can be any integer, and the factor should be independent of the 
choice of r. Although this fact is not manifest from Eq. (146), we can check it 
numerically, or analytically by expanding it in terms of e. 

We thus have two expressions for the ingoing- wave function i?™. One is given 
by Eq. (97) with pf^ expressed in terms of a series of hypergeometric functions 
as given by Eq. (101), which converges everywhere except at r = oo. The other 
is expressed in terms of a series of Coulomb wave functions given by 



R" 



(147) 



which converges at r > r+, including r = oo. Combining these two, we have a 
complete analytic solution for the ingoing-wave function. 

Now we can obtain analytic expressions for the asymptotic amplitudes of 
i?'", B^"'^'^^, B™^ and B'^^^ . By investigating the asymptotic behaviors of the 
solution at r — > oo and r — > r+, they are found to be 



_gref 



= W 



Jn 

) 

-J7r,.sin7r(;/ - s + ie) 



n— — OO 



-l-2s 



sin7r(i^ + s — iej 



— it 111 



(148) 

(149) 
(150) 
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Incidentally, since we have the upgoing solution in the outer region (140), 
it is straightforward to obtain the asymptotic outgoing amplitude at infinity 
(^trans f^^^ ^q. (134). Ws find 




(151) 



4.5 Low frequency expansion of hypergeometric expan- 



So far, we have considered exact solutions of the Teukolsky equation. Now let 
us consider their low frequency approximation and determine the value of v. 
We solve Eq. (118) with n = 0, 



with a requirement that v ^ £ as e ^ 0. 

To solve Eq. (152), we first note the following. Unless the value of v is such 
that the denominator in the expression of a^-^ or 7^ happens to vanish or (3'^ 
happens to vanish in the limit e ^ 0, we have a'^ = 0(e), 7^ = 0(e) and 

= 0(1). Also, from the asymptotic behavior of the minimal solution as 
n ±00 given by Eq. (116), we have Rn{v) = 0{e) and L^n{v) = 0(e) for 
sufficiently large n. Thus, except for exceptional cases mentioned above, the 
order of a!^ in e increases as |n| increases. That is, the series solution naturally 
gives the post-Minkowski expansion. 

First let us consider the case of Rn{v) for n > 0. It is easily seen that 

= 0(e), 7;; = 0(e) and = 0(1) for all n > 0. Therefore, we have 
Rniv) = 0{e) for all n > 0. 

On the other hand, for n < 0, the order of L-n{v) behaves irregularly for 
certain values of n. For the moment, let us assume that L-i{v) = 0(e). We see 
from Eqs. (106) that afj = 0(g), 7^ = 0(e), since v = £ + 0(e). Then, Eq. (152) 
implies (}q = O(e^). Using the expansion of A given by Eq. (92), we then find 
u = (.-{- O(e^), i.e., there is no term of 0(e) in v. With this estimate of v, we 
sec from Eq. (110) that L-i{v) = 0(e) is justified if L_2(z^) is of order unity or 
smaller. 

The general behavior of the order of L_„(z/) in e for general values of s is 

rather complicated. However, if wc assume s to be a non- integer and (. > |,s|, 
and T = {e — mq)/K, = 0(1), it is relatively easily studied. With the assumption 
that V = i + O(e^), we find there are three exceptional cases: For n = —21 — 1, 
we have Q!„ = 0(e), /3 = O(e^) and 7„ = 0(e). For n = — £ — 1, we have 
aj; = 0(l/e) and = 0(l/e) and 7„ = 0(e). And for n = we have 
a„ = 0(e), j3n ~ 0(l/e) and 7„ = 0(l/e). These implies consecutively that 
L_2e-i{i') = 0(l/e), L_£_i{u) = 0(1) and L_^{u) = O(e^). To summarize, we 
have 



sion 



A';+a;^i?i+7oi-i =0, 



(152) 



J n 



= 0(e) for all n > 



fn-l 



35 



I/_2^-l(z^) 
Ln{v) 



f-e 



f-e+1 

/-2€-l 



= 0{e% L_,_i(^) = %i = 0(l), 



f-e 



n 



0(l/e), 



2t 



J n 



J n- 



0{e) for all the other n < 0. 



(153) 



n+l 



With these results, we can calculate the value of u to 0(e), which is given by 



1 



21+1 



+ 



[(^ + 1)^ 



„2l2 



{2£ ~ l)2£{2i + 1) 



£{£ + !) (2£+ 1)(2£ + 2)(2£ + 3) 



(154) 



Now one can take the limit of an integer value of s. In particular, the above 
holds also for s = 0. Interestingly, v is found to be independent of the azimuthal 
eigenvalue m to O(e^). 

The post-Minkowski expansion of homogeneous Teukolsky functions can be 
obtained with arbitrary accuracy by solving Eq. (105) to a desired order and 
by summing up the terms to a sufficiently large \n\. First few terms of the 
coefficients are explicitly given in [43]. A calculation up to a much higher 
order in 0(e) was performed in [71], in which the black hole absorption of 
gravitational waves was calculated to 0{v^) beyond the lowest order. 



5 Gravitational waves from a particle orbiting a 
black hole 

Based on the ingoing wave functions discussed in the previous two sections, we 
can derive gravitational wave energy and angular momentum flux emitted to 
infinity. The formula for the energy and the angular momentum luminosity 
to infinity are given by Eqs. (35) and (36). Since most of the calculations are 
very long, we only show the final results. In [46], some details of the calcula- 
tions are summarized. We define the post-Newtonian expansion parameter by 
X = (Mfiip)-'^/^ where M is the mass of the black hole, is the orbital angu- 
lar frequency of the particle. Since the parameter x is directly related to the 
observable frequency, this result can be compared with the results by another 
method easily. 



5.1 Circular orbit around a Schwarzschild black hole 

First, we show gravitational waves luminosity emitted by a particle in circular 
orbit around the Schwarzschild black hole [73, 77] . In this case, Q^, is given, in 
the Schwarzschild coordinate, by fl^ = (M/rg)^/^ = Oc where ro is the orbital 
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radius. 
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where {dE/dt)^ is the Newtonian quadrupolc luminosity given by 
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(156) 



This is the 5.5PN formula beyond the lowest, Newtonian quadrapole formula. 
We can find that our result agrees with the standard post-Newtonian results up 
to 0{x^) [14] in the limit /i/M < 1. 

5.2 Circular orbit on the equatorial plane ciround a Kerr 
black hole 

Next, we consider a particle in circular orbit on the equatorial plane around a 
Kerr black hole [74]. In this case, the orbital angular frequency is given by 



(157) 



where ilc is the orbital angular frequency of the circular orbit in Schwarzschild 
case, V = (M/ro)^^^, q ~ a/M, and ro is the orbital radius in the Boycr- 
Lindquist coordinate. The effect of the angular momentum of the black hole is 
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given by the corrections depending on the parameter q. Here, q is arbitrary as 
long as |g| < 1. The luminosity is given up to O(x^) (4PN order) by 
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5.3 Slightly eccentric orbit around a Schwarzschild black 
hole 

Next, we consider a particle in slightly eccentric orbit on the equatorial plane 
around a Schwarzschild black hole ([46], Sec. 7). We define ro as the minimum of 
radial potential R(r)/r'^. We also define an eccentricity parameter e from the 
maximum radius of the orbit Tmax, which is given by fmax = ^^{l + e). These 
conditions are explicitly given by 



d{R/r^) 



0, and i?(ro(l + e)) =0. (159) 



dr 

We assume e <C 1. In this case, is given to 0{e'^) by 



r=ro 

2\ 



where Q,c = (M/ro)^/^ is the orbital angular frequency in the circular orbit case. 

We show the energy and angular momentum luminosity which are accurate to 
0{e^) and to 0{x^) beyond the Newtonian order. They are given by 
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where {dJ/dt)^ is the Newtonian angular momentum flux expressed in terms of 

X, 



~dt 



32 f fi 
M 
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and the e-independent terms in both (dE/dt) and {dJ/dt) are the same and are 
given by the terms in the case of circular orbit, Eq. (155). 



5.4 Slightly eccentric orbit around a Kerr black hole 

Next, we consider a particle in slightly eccentric orbit on the equatorial plane 
around a Kerr black hole [70]. We define the orbital radius Tq and the eccen- 
tricity in the same way as in the Schwarzschild case by 



d{R/r'^ 



dr 



= 0, and R{ro{l + e))=Q. 



(164) 



We also assume e <C 1. In this case, Q,^ is given to O(e^) by 
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We show the energy and angular momentum luminosity which are accurate 
to O(e^) and to O(x^) beyond Newtonian order. 
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5.5 Circular orbit with small inclination from the equato- 
rial plane eiround a Kerr black hole 

Next, we consider a particle in circular orbit with small inclination from the 

cqiiatorial plane around a Kerr black hole [69]. In this case, apart from the 
energy E and z-component of the angular momentum 1^, the particle motion 
has another constant of motion; the Carter constant C. The orbital plane of 
the particle processes around the symmetric axis of black hole, and the degree 
of precession is determined by the value of the Carter constant. We introduce 
a dimensionless parameter y defined by 



y=^' Q''=ll+a\l-E% (168) 



Given the Carter constant and orbital radius tq, the energy and angular mo- 
mentum is uniquely determined by R{r) = 0, and dR{r)/dr = 0. By solving 
the geodesic equation with the assumption y <C 1, we find that y^/^ is equal to 
the inclination angle from the equatorial plane. The angular frequency fl^ is 
determined to 0{y) and O(u^) as 



(169) 

Wo show the energy and the ^;-component angular flux to 0{v^). 
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Using Eq. (169), we can express v in terms of a; = (Mfi^)^/^ as 
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We then express Eqs. (170) and (171) in terms of x as 
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5.6 Absorption of gravitational waves by a black hole 



(174) 



In this section, we evaluate the energy absorption rate by a black hole. The 
energy flux formula is given by [79] 
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In calculating Zf^^, we need to evaluate the upgoing solution i?"P, and 
asymptotic amplitude of ingoing, and upgoing solution, B'""^, B^'^^'^^, and C'*''*"^ 
in Eqs. (14) and (15). Evaluation of the incident amplitude iJ*''^"'' of the ingoing 
solution is essential in the calculation. Poisson and Sasaki [59] evaluated them, 
in the case of circular orbit around the Schwarzschild black hole, up to 0(e) 
beyond the lowest order, and obtained the energy flux at the lowest order, 
using the method we have described in section 3. Later, Tagoshi, Mano, and 
Takasugi [71] evaluated the energy absorption rate, in the Kerr case, to 0{v^) 
beyond the lowest order using the method in section 4. Since the resulting 
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formula is very long and complicated, we only show here to 0{v^) beyond the 
lowest order. The energy absorption rate is given by 
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and ■(/''"H^) is the polygamma function. We see that absorption effect begins at 
0(v^) beyond the quadrapole formula in the case q 0- If we set g = in the 
above formula, we have 



(179) 



which was obtained by Poisson and Sasaki [59]. 

We note that the leading terms in {dE/dt)n are negative for q > 0, i.e., the 
black hole loses the energy if the particle is co-rotating. This is because of the 
superradiance for modes with A; < 0. 



6 Conclusion 

In this article, we described analytical approaches to calculate gravitational 
radiation from a particle of mass /x orbiting a black hole of mass M with M ji, 
based upon the perturbation formalism developed by Teukolsky. A review on the 
this formalism was given in section 2. The Teukolsky equation, which governs 
the gravitational perturbation of a black hole, is too complicated to be solved 
analjrtically. Therefore, one has to adopt a certain approximation scheme. The 
scheme we employed is the post-Minkowski expansion, in which all the quantities 
are expanded in terms of a parameter e = 2Muj where uj is the Fourier frequency 
of the gravitational waves. For the source term given by a particle in bound 
orbit, this naturally gives the post-Newtonian expansion. 

In section 3, we considered the case of a Schwarzschild background. For 
a Schwarzschild black hole, one can transform the Teukolsky equation to the 
Regge- Wheeler equation. The advantage of the Regge- Wheeler equation is that 
it reduces to the standard Klein-Gordon equation in the flat space limit, and 
hence it is easier to understand the post-Minkowskian or post-Newtonian eifects. 
Therefore, we adopted this method in the case of a Schwarzschild background. 
However, the post-Minkowski expansion of the Regge- Wheeler equation is not 
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quite systematic, and as one goes to higher orders, the equations to be solved be- 
come increasingly complicated. Furthermore, for a Kerr background, although 
one can perform a transformation similar to the Chandrasekhar transformation, 
it can be done only at the expense of losing the reality of the equation. Thus, 
the resulting equation is not quite suited to analytical treatments. 

In section 4, we described a different method, developed by Mano, Suzuki 
and Takasugi [45, 43], that directly deals with the Tcukolsky equation and con- 
sidered the case of a Kerr background with this method. Although the method 
is mathematically rather complicated and it is hard to obtain physical insights 
into relativistic effects, it has a great advantage that it allows a systematic 
post-Minkowski expansion of the Teukolsky equation, even on the Kerr back- 
ground. We gave a thorough review on how this method works and how it gives 
a systematic post-Minkowski expansion. 

Finally, in section 5, we recapitulated the results of calculations of the grav- 
itational waves for various orbits that had been obtained by various authors by 
the methods described in sections 3 and 4. These results arc useful not only by 
themselves for the actual case of a compact star orbiting a supermassive black 
hole, but also give us useful insights into higher order post-Newtonian effects 
even for a system of equal-mass binaries. 
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